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Abstract 

We present a detailed computation of the fully exclusive cross section of p + p ^ + 
7 + X with X = and 1 jet in the framework of the factorization theorem and dimensional 
regularization. Order as and photon bremsstrahlung contributions are discussed in the 
MS mass factorization scheme. The resulting expressions are ready to be implemented 
numerically using Monte Carlo techniques to compute single and double differential cross 
sections and correlations between outgoing pairs of particles. 
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I INTRODUCTION. 

Ever since Mikaelian's discovery of a zero in the amplitude of the partonic subprocess q + q —>■ 
W+'-f , radiative production of W bosons has been discussed as a way of testing the vahdity 
of the Electroweak Theory. The study of differential distributions in p + p ^ W + 'y + X 
may be the best way to place bounds on the magnitude of the magnetic dipole and electric 
quadrupole moments of the W boson. Deviations from the Standard Model could show up as 
a shift of the photon distributions near the dip that is a reflection of the partonic zero. The 
QCD corrections to the reaction p + p—^W + 'j + X and its deviations from the Standard 
Model have been studied in and other references therein. These papers have 

been mainly devoted to the analysis of single photon distributions, photon-W boson pair 
mass correlations and charged lepton-photon pseudorapidity correlations and they either 
neglect or approximate the photon bremsstrahlung contributions. 

When computing the photon inclusive process in 0] and all the singularities associated 
with a jet emitted in a coUinear or a soft region of phase space were regularized by analyt- 
ically performing all the integrals associated with the jet and the W boson in n space-time 
dimensions. Although the numerical advantage of this procedure is obvious -one is left with 
only photon variables to be integrated over numerically- the predictive power of the whole 
computation is limited by the fact that one looses information about the energies and angles 
of the jet and the W boson. 

In the present work we redo the exact first order calculation reported in in an exclu- 
sive fashion. We present analytical results for the integrands needed in the computation of 
physical observables related to any of the outgoing particles in the reactions p + p — > W^"*" + 7 
and p + p — > + 7 -|- jet. Using these results we will extend the studies of the Elec- 
troweak and QCD sectors of the Standard Model by providing a complete set of single and 
double differential distributions and correlations including the W boson and, when applica- 
ble, the jet. Deviations of the experimental data from the theoretical predictions could not 
only mean new physics in the Electroweak sector, but would also probe the QCD behavior 
and the underlying photon bremsstrahlung processes. In particular, an inadequate photon 
bremsstrahlung approximation would also result in deviations from the predicted photon 
single and double differential distributions and correlations. 

The method that we employ for computing exclusive cross sections is based on the one 
used by Mele et. al. in the context of pair production and production of heavy quarks 0. 
This method allows for control of all soft and initial (final) state coUinear singular regions of 
phase space in the framework of dimensional regularization and the factorization theorem. 

We consider three different scenarios: (1) the 2-body inclusive production of and 7, 
(2) the exclusive production of W~^, 7 and 1 jet and (3) the exclusive production of and 7 
accompanied by jet. In all three cases we take into account exact 0{as) QCD contributions 
in the MS mass factorization scheme. Contributions arising from photon-quark and photon- 
gluon fragmentation functions (generically called "photon bremsstrahlung contributions") 
are also included in our discussion. 

In Section II we present a detailed review of the process p+p — > + '-f + X for X = 0, 1 
jet in the framework of the parton model and the factorization theorem. 

In Section III we show how the cancellation of singularities is performed in an exclu- 
sive fashion in each of the hard scattering channels of our process in the framework of n 
dimensional regularization. 
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Section IV is devoted to the definition of tlie tliree experimental scenarios and tlieir 
corresponding cuts. 

We end our study in Section V witli a discussion of tlie numerical implementation of the 
several expressions for the cross sections, together with a listing of the relevant formulae. 

Results for total, single and double differential cross sections and correlations between 
pairs of outgoing particles are given in a separate paper . 
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II THE PROCESS p + p ^ W+ ^ +X, THE PAR- 
TON MODEL AND THE FACTORIZATION THE- 
OREM. 

II.A Introduction. 

We are considering the hadronic processes given by 

p + p-^Vr+ + 7, (2.1) 

and 

p + p + ^ + jet. (2.2) 

In what follows we will omit the charge index when referring to the boson. In 
the framework of the parton model we can formally write the hadronic 2-body inclusive 
differential cross section in the following way: 

^ nn^'nn [^(^^^ ^"^^^^ ^^^^^ ^^^^^ ^] = 

V f / dui [ du2 f du3 Dip{ui)Dyp{u2)D.,k{u3) 
TTi. \Jo Jo Jo 



E 



- DQ^DQ^ 



(2.3) 



where J2x is a sum over sets X of physical particles integrated over their phase space. 
J2i,j,k denotes sums over partons k (by partons we mean quarks, antiquarks, gluons and 
photons.) J2x is a sum over sets x of outgoing partons integrated over their phase space. In 
n dimensional space-time DQi and DQ2 are generically given by: 

where Qi is an n-momentum vector with space components Qi and time component Qi^. 
Y^^D^a^ / DQ1DQ2 formally denotes the "bare" partonic 2-body inclusive differential cross 
section, which can be directly computed using perturbation theory in the Standard Model. 
DiAiu) for A = p,p are the bare partonic densities. The bare fragmentation function D^kiu) 
gives the photon momentum fraction density when a parton of type k and momentum q 
fragments into a photon of momentum uq and any number of hadrons. After renormaliza- 
tion has been performed the bare partonic cross sections in (2.3) contain soft and collinear 
singularities coming from virtual corrections as well as integration over phase space of non- 
empty sets of outgoing partons x. Cancellation of soft singularities will occur after addition 
of the soft pole terms in the virtual corrections with the soft pole terms in the corresponding 
emission processes. By virtue of the factorization theorem ||^ the bare partonic densities and 
the bare fragmentation function are defined to contain singularities that cancel against the 
remaining collinear singularities in the bare partonic cross sections so that the hadronic cross 
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section on the LHS of (2.3) is a finite quantity. This procedure is implemented at a specific 
mass factorization scale M . We will define this scale to equal the renormalization scale /i. To 
avoid complicating the subsequent formulae this scale dependence is not explicitly written, 
except where necessary for the discussion. According to the factorization theorem we can 
rewrite the singular bare partonic 2-body inclusive cross section in terms of non-singular 
"hard scattering cross sections": 

"1 rl rl 



E {I dVll "^""'i <i (^l)<i(^2)t^fee(^3) 

Q2, 



12 



\ DQ,DQ2 



a{vipi) b{v2P2) W{Qi) c(- 

^^3 



(2.5) 



where J2y D^(t / DQ1DQ2 denotes a 2-body inclusive hard scattering differential cross section. 

d^i{v) denotes the splitting function of parton i into a parton a and a set of par- 
tons Xa with a carrying a momentum fraction v of its parent parton i. These splitting 
functions factorize the coUinear singularities contained in the bare partonic cross section 
Y^r^ D^a^ / DQ1DQ2 and they can be exactly computed order by order in perturbation 
theory. In this way (2.5) is solved perturbatively for the hard scattering cross sections 
j:yD^a/DQ,DQ2. 

Using (2.5) in (2.3) we can rewrite the hadronic 2-body inclusive differential cross section 
in terms of only non-singular quantities: 

E nn^'nn t^^^^) P^^^^ ^ ^^^^^ ^^^2) ^] = 

"1 rl rl 



nil / '^^2 / dT'i fapiTi)fbp{T2)fjciT3) 

"T, \J0 Jo Jo 



X 



E 



D 



a 



a{nP,) b{T2P2) ^ W{Q,) c(^) X 

7-3 



(2.6) 



DQiDQ2 

The parton densities fap^T), fbp{T) and the fragmentation function /7c (t) are defined by: 

'■1 1 ,T. 



fap{r) = I du -dlt{-)Dip{u) 

hp{r) = tdu-dt^{-)D 
^ ■ Jt u u 



jp\U) 



Mr) = E f du-D,,{u)dll{-), (2.7) 

XcfK 

and are obtained by fitting data of deep inelastic scattering to results of perturbation theory 
at a given mass factorization scale M. At present there is not enough data available to fit 
the photon fragmentation functions so one has to rely upon an approximation, for example, 
the so called leading- log-approximation P|,[pi]|]. 

In the computation of hadronic quantities we use (2.6) with 
a, 6 e {g, g, g} and c e {g, g, 7} where g, g and g denote quark, antiquark and gluon respec- 
tively. The photon (7) is treated in a dual way: it is a hadron, i.e. an observable final state 
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particle, and it is also a parton of our Lagrangian. In our 1-body inclusive computation in 
we only considered contributions from /^^ and neglected f^c for c e{q,q,g}. In the present 
work we include the four contributions keeping terms up to 0{asaaw), where as, a and aw 
are the strong, electromagnetic and electroweak fine structure constants. 



II.B Contributions from /^^. 

The leading photon-photon splitting function is given by dj^^{T/u) = 6(1 — r/u), i.e. when 
no partons are emitted from the photon. This leading order splitting function can be iden- 
tified with the bare fragmentation function D^^{u) when zero hadrons are fragmented from 
the photon. Using this in (2.7) we obtain the leading contribution to the photon-photon 
fragmentation function 

/,,(r) = 5(l-r). (2.8) 
Setting c = 7 in (2.6) and keeping hard scattering contributions up to 0{asaaw) we obtain 

dn £ dT2 lu{T,)Mr2) [ j^Q^^Q W^Pi) <li^2P2) - W{Q,) 7(Q2)] 



a 



DQiDQ, 



[g(riPi) g(r2P2) ^ W{Q{) 7(Q2) g] 



+/,p(ri)/5,(r2) -^^^^ WiPi) q{r2P2) - W{Q,) ^{Q^) g] 

+ {q ' — ^ q) } • 

(2.9) 

Note that the lowest order hard scattering cross section is always equal to the corresponding 
lowest order bare partonic cross section, as we will verify in subsections F,G and H. 



II. C Contributions from f^q and f^q. 

Setting c = g and c = g in (2.6) we obtain for these contributions: 



E 

X 



^2„H \ 79 



DQ1DQ2 



b(Pi) p(P2) ^ W{Q,) 7(Q2) X] 



dn / dT2 / dr- 



X <^ fqp{Ti)fgp{T2)f^g{T3) 



+ fgp{.ri)fgp{T2)f^q{n) 



DQ^DQ2 
DQ^DQ2 



q{nP,) g{T2P2) ^ I^(Qi) g( 



Q2 

1"3 



g{nP,) g(r2P2) ^ ly(Qi) g(^) 



(2.10) 



and 



-a 



^ \DQ1DQ2 



b(Pi) piP2) ^ W{Q,) 7(g2) X] 



JO JO 



X ygpin)fqpir2)f^g{T3) 
+ fqp{Tl)fgp{T2)f^g{T3) 



DQiDQ2 
DQ^DQ2 



,Q2 



g(riPi) (7(r2P2) - W{Q,) g(^) 



(2.11) 



II.D Contributions from f^g 

Setting c = g m (2.6) we have: 



^ \DQ1DQ2, 



[p(Pi) p(P2) ^ w^(Qi) 7(^2) X] = [' dn t dr2 C dn 

Jo Jo Jo 



X yqpiTi)fgp{T2)f^g{T3) 



DQiDQ2 



qinP^) q{r2P2) ^ W{Qi) g 



Q2 



(2.12) 



In B, C and D sums over flavors of quarks q and antiquarks q satisfying tlie electric 
cliarge conservation for production of are implicit. For a more detailed discussion on 
the way we treat this issue see section IV of [Q. The contributions from C and D will be 
referred to as "photon bremsstrahlung contributions" . 



II. E The incoming hard scattering channels. 

We can now regroup all terms in B, C and D according to three partonic channels in the 
incoming hard scattering state: 



X i^U{n)fqp{T2) {^ j^qI^q^ WiPi) q{r2P2) - W{Q,) 7(g2)] 



+/o ^-^"^^DgT^g, 



q{TiPi) q{T2P2) ^ W{Q,) g{^] 

T3 



(2.13) 



+ [ dri [ dT2 I dr-^ 
Jo Jo Jo 

X { fqp{n)fgp{T2)f^q{rz) 



3 



DQiDQ2 



73 J 



(2.14) 



X |/,p(ri)/,,(r2) b(nPi) g(r2P2) ^ l^(gi) 7(^2) 

+/,p(ri)/,p(r2) j^qI^q^ WiPi) g^r^P,) ^ W^(gO ^{Q,) q] 



+ [ dri [ dT2 [ dr^ 
Jo Jo Jo 

X <^ fgp{n)fqp{T2)f^q{n) 



3 



DQ^DQ^ 



DQ1DQ2 



g{TiPi) q{T2P2) ^ W{Qi) g(^) 

T3 . 

g(TlPl) ^(T2P2) ^ W{Q,) q{^) 

73 . 



(2.15) 



II.F Factorization in the channel. 

Let us write (2.5) for i — q, j = q and k = j: 

E ^^^^ fe) - w^(Qi) 7(Q2) x] = 



E dv2 £ dv3 dll{v^)dll{v2)d^^l{v3) 
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a(t^iPi) h{v2P2) W{Qi) c(— ) ?/ 



(2.16) 



We will solve ( |2.16| ) at 0{aaw) and 0{as(y.aw) so only x = { } and x = {g} contribute on 
the LHS. By constraining the sums on the RHS so that x^ U U Xc U ?/ C a; we obtain at 
0{a%): 



[q{Vi) q{P2) ^ W{Q,) ^{Q^) 



DQ1DQ2 



[qiPi) q{P2) ^ W{Q{) 7(Q2)] 



(2.17) 



and at 0{as)'- 



DQiDQ2 



[q{pi) q{P2) ^ W{Q,) 7(Q2)] 



DQ^DQ, 



[q{pi) q{P2) ^ W{Q{) 7(g2) 9] 



[q{pi) q{P2) ^ W(Qi) im] 



DQ1DQ2 

+ \q{Pi) q{P2) - Wm 7(g2) g] 



DQ1DQ2 

'dv PmJ^^^ [q{:vpi) q{p2) ^ W{Q,) 7(Q2)] 



as 



27re Jo 



"5 



27re Jo 



"''^'DQiDQ2 
1 _ dV(°) 



[g(Pi) g(^^P2) ^ H^(gi) 7(^2)] 



(2.18) 



and the corresponding equations for pi < — > p2. In deriving ( p.l7|) and ( p.l8|) we have used 
the following splitting functions: 



dll\v) 



"bq 



27re 

as 

27re 



-Pqq{v)5aq 
Pqq{v)6bg 



(2.19) 



and the definition: 



with 



3 = i - 7s + ln(47r) 
e e 



(2.20) 



(2.21) 



and e = (4 — n)/2. The running strong fine structure constant is as = 5'|(/i)/47r. In the 
MS mass factorization scheme Kij{v) = for all relevant i,j in this and the following two 
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subsections. 



For i = j = q{q) we have: 



Cf 
Cf 



- V 



with 



'dv\-^] f{v) 



dv + dv 

li — V ) Jo 1 — V Jvo 1 — V 

where < f o < 1 and the color factor is given by Cf = 4/3. 

In an analogous way, setting i = q, j = q, k = g in (|2.5|) we obtain: 



(2.22) 
(2.23) 



[qiPi) q{P2) ^ W{Qr) giQ^)] 



DQ^DQ2 



[q{pi) q{P2) ^ W{Q{) g{Q2)] . 



(2.24) 



II. G Factorization in the qg channel. 



Setting i = q,j = g and A; = g in ( p.5|) and again keeping terms up to 0{asaaw) we can 
only have x = { }, thus constraining the sums on the RHS to U a;;, U U ?/ = { }. We 
obtain at 0{as)- 



DQ1DQ2 



[q{pi) 9{P2) - W{Q{) q{Q2)] 



DQ1DQ2 



[q{pi) 9{P2) -> W{Q,) q{Q2)] . 



(2.25) 



Resetting A; = 7 in (|2.5| ) we can now only have x = {q}, thus constraining the sums on the 
RHS to Xa l-i Xb U Xc U y = {q}. We obtain at 0(0:5): 



DQ1DQ2 



[q{pi) 9{P2) ^ W{Q,) 7(g2) q] 



DQ1DQ2 



[q{pi) 9{P2) ^ W{Qr) ^{Q^) q] 



as 
2Txe Jo 



dv Paniv 



, dv P^niv) 

27re7o ' DQ1DQ2 



a 



[q{pi) q{vp2) ^ W{Qr) ^{Q^)] 
q{pi) 9{P2) - W{Q,) g(^ 



(2.26) 

and the corresponding equations for pi < — > p2- The following splitting functions have been 
used when deriving ( p.26|) : 



27re 

Oi 



(2.27) 
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with 



+ {1- vf 



2i + {i-vy 



(2.28) 



where Cg = — 1/3 is the charge of the outgoing quark q on the LHS of ( p.26| ), in units of e 
and a = e^(/i)/47r is the running electromagnetic fine structure constant. 



II.H Factorization in the gq channel. 

Analogously to the previous case, setting i = g, j = q and /c = g, 7 in ( p.5|) and keeping up 
to 0{as<y<yw) we obtain: 



DQ^DQ, 



DQ^DQ2 



(2.29) 



and 



DQiDQ2 



Wi) q{P2) ^ W{Q,) 7(g2) q] 



[giPi) q{P2) ^ W{Q{) 7(g2) q] 



as 



a 



27re Jo 



[q{vp,) q{p2) ^ ly(Qi) 7(^2)] 
giPi) q{P2) - W{Q{) g(^) 

V 



(2.30) 



and the corresponding equations for pi < — > p2. In deriving ( ^.3(J| ) we have used the following 
splitting functions: 



with 



a 



^2 + (1 - vf 



(2.31) 



l + (l-t;) 



(2.32) 



where Cg = —2/3 is the charge fraction of the outgoing antiquark q on the LHS of ( p.30| ), in 
units of e. 
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III EXCLUSIVE CANCELLATION OF SINGULAR- 
ITIES. 



III.A Introduction. 



In order to compute the hard scattering cross sections required in ( p.l3| )-( pT5| ) using ( |2.17|) , 
QCTI ), p:^)-(P^), (!r2% and d^;^) we first need the bare partonic cross sections D'^a^ /DQ1DQ2 
on the RHS of these equations evaluated in n = 4 — 2e space-time dimensions. 

When computing phase space integrations of outgoing massless particles singularities 
appear in regions of phase space where one of these particles is collinear to any other massless 
on-shell parton or where one of the outgoing massless gauge bosons is soft. Since we are 
tagging the outgoing photon, we don't have to worry about singularities associated with 
integration over the photon's phase space. 

With this in mind we will classify the 2 to 3 body Feynman amplitudes according to the 
way the outgoing massless particle q,q or g which is integrated over is attached to the rest of 
the legs of the diagram. In Fig. 1 we have decomposed the 2 to 3 body Feynman amplitude 
for the partonic reaction q + q—>-W + '-f + g into three pieces: M'^^^'^'^^ = Mj^+Mj^+Mjjj: 



q(pi) 




q(Pi) 




q(Pi) 



q(Pi) 



q(P2) y 

j^qq^Wyg 





+ 



q(P2) y q(P2) y q(P2) 

Figure 1. Decomposition of the 2 to 3 body Feynman amplitude in the qq channel. 



The labels "g" near the solid vertices at the end of the incoming quark and antiquark legs 
in the amplitude Mjjj mean that these vertices do not contain the outgoing gluon. Arrows 
show the direction of the fermionic charge and the charge. The momenta pi and p2 are 
always incoming. The shaded blobs denote the inclusion of all possible Feynman diagrams 
(except for the mentioned constraints in Mjjj.) If we integrate the squared matrix element 

over the phase space of the outgoing gluon summed over physical polarizations the J2 



and pieces of the squared matrix element have collinear and soft singularities while 

the interference term [J2 Mj^Mj^'^ + c.c) has only soft singularities. Other pieces of the 
squared matrix element have no singularities. 

For the partonic reaction q + g^W + 'y + qwe have M''^'^^^^ = M]^ + Mff + Mff^, 
as shown in Fig. 2. 
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q(Pi) 




q(Pi) 



g(P2) 



jyjqg^Wyq 




q(Pi) 




q(Pi) 




g(P2) 



g(P2) 



Figure 2. Decomposition of the 2 to 3 body Feynman amplitude in the qg channeL 

In this case, only the E \M1^\' and E |M/f T pieces of the squared matrix element have 
coUinear singularities, with no soft singularities in this channel. 

In Fig. 3 we show the analogous decomposition for the partonic reaction g + q 
W + + q: Ms^"^^T5 = Mff + Mff + Mffj. As in the previous case, only the E Mff 



and E 



Mff 



pieces of the squared matrix element heve coUinear singularities and no soft 



singularities are present in this channel. 



g(Pi) 



g(Pi) 






g(Pi) 



q(P2) y 

.gq-^ Wyq 



q(P2) 



q(P2) y q(P2) y 

M^^^"^^^ M§ Mif xvxni 

Figure 3. Decomposition of the 2 to 3 body FejTiman amplitude in the gq channel. 




We will develop the above decompositions in more detail for each channel in subsections 
III.B, C and D. 

The 2 to (2 + /) body partonic differential cross section for the W + 'y + X production 
process is defined in n dimensions of phase space by: 



^^{2+1) 



1 1 



N2s 

xEi^r 



DQ^DQ^Dh ■ ■ ■ Dh{2'KY5^{pi + P2 - Qi - Q2 - h) 



(3.1: 



where we have averaged over possible incoming states of different polarizations and colors 
and the sum on the RHS is over polarizations and colors of all particles. Qi and Q2 are the 
n-momenta of the boson and the photon respectively. The variables in (|3.1| ) are defined 
in the following way: 



s = 2pi- P2 

^ n-2 
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DQ2 
Dki 



1 



2(27r)"-i 
Dk 



Q2 
1 



n—3 



d 



Q2 dVtn-2 (Q 



2 [.(Pl, • • • , (Pn-2 



2(27r)"-i 



n—3 



d^ln-2 [k (^1, . . . , V'n- 



(3.2) 



The angle differentials in ( p.2| ) are generically given by: 

(p • • • 5 ttn-2)) = rfcosai sin"'~'^ai ■ ■ -dcosan-a sin" a„_3(iQ;„_2 



(3.3) 



To account for the experimental cuts on the outgoing particles that define the experimental 



scenario under consideration we have to include on the RHS of ( |3.1| ) an extra factor of 
C{Qii Q2, ki, . . . , ki). These cuts may be expressed in a covariant way in terms of G step 
functions as we will see in more detail in section IV. In the rest of this chapter and the 
following sections we omit the charge index when referring to the boson. 

III.B The qq CHANNEL. 

The 2 to 2 body total partonic cross section for the reaction qq —>■ W'j may be written in 
the following way: 



DQ1DQ2 



2^P 



[q{pi)q{p2) ^ W{QMQ2)] 



DQ^DQ2 



4iV2 



(3.4) 



with 



$(.) 
f3{s) 
Pis) 



)2e 



r(i-e) 

l-p(.) 



Any 1 

s J IQtc 



l-2e/ 



The only independent invariants in the squared matrix elements in 

b = 2p,-Q2 = ^Pis){l + cosei). 

The rest of the invariants may be expressed in terms of b and s as follows 

s-b 



are s and 



(3.5) 



(3.6) 



2pi ■ Qi 
2p2 ■ Qi 
2p2 ■ Q2 
2Qi ■ Q2 



b + Mlr 

s~M^-b 

s-M^. 



(3.7) 
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In ( p.4|) we chose the {n — l)-th axis to be the one pointing in the direction of pi. The 
integration over the angles 6*2, ... , 9n-2 has been performed because there is no dependence 
on them in the 2 to 2 body squared matrix element. To account for the experimental cuts 
on the outgoing particles we implicitly include an extra factor of C{Qi,Q2,0) on the RHS 
of ( |3.4|) . In ( |3.5|) and subsequent equations the unprimed variables refer to variables in the 
center of mass system of the incoming partons. 

To obtain the 2 to 3 body partonic cross section for the reaction qq — *■ W'jg we define a 
primed reference frame in the W'j center of mass system such that the n-momenta are given 
by: 



k' 

Q2 



(1,0,..., 0,0,1) 

P2,o(l>0,...,0,sinr/',cosr/' 
/cq (1, 0, ... , 0, sinip', cosip') 

Qi,Q 




. . . , sin 6[ sin 62 cos 6*3, sin 6[ cos 62, cos 6[ 
1, . . . , — sin 6^^ sin 62 cos 9'^, — sin 9[ cos 62, — cos 9[] 



Q'l 



{3.t 



The 2 to 3 body total partonic cross section in the qq channel may thus be written as follows: 



DQ,DQ2 



a 



DQiDQ2 
1 1 



Am 2s 

X 



dx ^{sx){l - xf-"^' dy(l-y^y 



(4vr) 



r(l-2e) 271 Jp{s) 
J^^d cos 9[ sm"^' 9[ £ d9'2 sin'^^ 0^ ^ |M'?^"^^t5(s, a, 6, c, rf) 



(3.9) 



where 



y = cos 

ko = (1 - x)pifi = (1 - a;) 



(3.10) 



We implicitly include on the RHS of ( p.9|) the factor C{Qi,Q2,k). We have chosen as 
independent invariants s and: 

a = 2p,-k='^{l-x){l-y) 

b = 2p^-Q2 = 2p[^^\Q[\{l + cos9[) 

c = 2k-Q2 = 2^0 Q[ (1 + sin ijj' sin 9[ cos 6*2 + cos ^p' cos 9[) 
d = 2p2 ■ Q2 = 2p2 Q'l (1 + sin rj' sin 9[ cos 6*2 + cos rj' cos 9[) 

(3.11) 
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so we have the following dependent invariants: 



2pi ■ Qi 
2p2 ■ Qi 

2gi ■ Q2 

2Qi ■ k 
2p2-k 



s — a — 

Mlr + a + h 
b + d ~ c 



= sx-M^ 
b-d 



s- M^ + c-a-b-d = -(1 



s 
2' 



x){l + y) 



(3.12) 



All invariants can now be expressed in terms of s, x, y, cos9[ and cos^^2 by solving for all the 
necessary quantities in the primed reference frame: 



Pi,o 

P2,0 

Qi,o 

COST]' 

cosip' 



/l + x+(l- x)y' 



X 



a/s /l + X — (1 — x)y' 



sx 



2 
Q'l 



1 + 



X 

sx 



sx 



P{sx) 



ix 



k' - ^ 



V(l+a;)2-(l 
1 — X + y{l + x) 
l + x + y{l 

^/s fl — x' 

2 



X 



X 



(3.13) 



Looking at (|3.9| ) it is clear that the soft divergences will be present in those pieces of 
squared matrix element that contain a factor (1 — x)~^, while the coUinear divergences will 
be due to factors of (1 — y)"^ or {1 + y)~^ in the squared matrix element. The squared 
matrix element for the qq channel can be written as: 



^ Im?^"-^^^'!' = ^ |Mff p + |Mf p + (J2 MZMjl* + c.c.) + remain . 



The type la matrix element (see Fig.l) can be written in the following way: 
Mia [g(pi, /i, Ai)g(p2, k) W{Q^)^{Q2)g{k, c, A)] = 



(3.14) 



2pi ■ k 



e'Jk) M [q {p, - k, l[, a) q{p2, h) - 1^(^1)7(^2)] 



(3.15) 



where li,l2 and l[ are quark color indices, Ai is the quark polarization index, c and A are 
the outgoing gluon color and polarization indices, a and p are Lorentz indices and gs is the 
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renormalized strong coupling constant. Other indices have been omitted because they are 
not necessary in what follows. The partial squared matrix element is given by: 



9l 



{2pi ■ kf 



xM* [q {p, - k, a') q{p2) ^ W{Q,)^{Q2)] . 



The Ih partial squared matrix element may be written in a similar fashion: 



Ek: 



^^^l\y CpR'^l,,Y.M[q{pMP2 -k,P)^ WiQMQ^)] 
xM* [q{pi)q{p2 - k,(3') W{QMQ2)] 



7130' 



(3.16) 



(3.17) 



The repeated indices on the RHS of ( |3.16| ) and ( |3.17| ) are contracted and the sum is over 
quark(antiquark) colors and polarizations as well as W and 7 polarizations. In the qq center 
of mass frame the tensors R^^'ja and R^^n, can be written in the following way: 



pact 

mJa 



R 



qqjb 



— 4pi ■ k 
-4p2 ■ k 



n 1 + y 
— + ^ 



7o 



2(1 -x) ; 

'2 — n ^ 1 — y 



2(1 -x) 



W+(i 



l + y 



1 — X 

1 — X 



^^1 



2(1 -x^ 
2(1 -x)^ 



/3'/3 



(3.18) 



In deriving ( |3.18|) we have summed over physical gluon polarizations in the covariant gauge 
where we may write: 



11-2 



A=l ■ 



(3.19) 



where, if = (fco, ^) then k = (^kQ,—lPj. The factors pi ■ k and p2 ■ k in front of the 
RHS of ( |3.18|) will cancel similar factors in the denominators of the RHS of (p.l6|) and 
( |3.17| ) respectively, leaving the type la and lb squared matrix elements with singular terms 
proportional to (1 — x)~^(l — y)""*^, {l~x)~^{l — y)~^ and (1 — x)~^(l +?/)~^ and (1 — x)~"'^(l + 
y)~^ respectively, i.e. both terms will contribute to the soft and collinear singularities. In a 
similar way the remain in ( |3.14| ) can be shown to have no collinear or soft singularities while 
the interference term M/aM/^ has only a soft singularity. It is thus convenient to define the 
non-singular function F'^^: 



F^^s,x,y,cose[,e',: 



4{pi-k){p2-k)J2\M'"'^^^' 



(3.20) 



so that we can now rewrite the 2 to 3 body total partonic cross section in (|3.9|) as follows: 



[qiPi)q{P2) W-fg] 



1 1 ,. ,,_2 r(i-e) s-'-' 



{Att 



r(l-2e) 271 
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X 



X 



j^^^ dx ^sx){l - x)-^-"^' ^dy (l - y"^) ^ ' J^^d cos 9[ sin-^' e[ 



de'^ sin-'' e'^ F"^ [s, X, y, cos e[, e'^) 



(3.21) 



Now that the singular factors in x and y have been isolated we can rewrite them as distri- 
butions for e < 0: 



'i-xY'-^' ~ 



{l-xo)-^'5{l-x) + 



2e 
(2yo 



1 — X 

1 



-2e 



f ln(l - 



X] 



2e 



[5{l-y)+5{l+y)] + 



\ 1 — X 
1 f 



2 \ l-y 



+ 



2 \ l + y 



+ 0(6) 



S/O 



(3.22) 



with C{2) = 7r^/6. We have introduced the following definitions: 

• 1 



1 

-1 



dxf{x) 



dyf{y) 



dyf{y) 



1 


1 


— X 




1 


1 


-y 




1 



+ 



1 + 2/ 



yo 



p(s) 1 — X 

dy{^ + 
■1 1 - y 



dy 



f{y) 



l+j/o 1 + 2/ 



I- X 

1-2/ 

-^+^»^^ /(2/)-/(-l) 
1 1 + 2/ 



(3.23) 



The parameters Xq and are arbitrary as long as they satisfy the conditions p{s) < Xq <1 
and < ?/o < 2. The symbol ~ in (|3.22|) means that the equality only holds under an 
integration over x ranging from p(s) to 1 for the first expression in ( p.22| ) and under an 
integration over y ranging from —1 to 1 for the second expression. When x and y are not 
integrated over their whole range care has to be taken when defining xq and y^ so as not to 
introduce unphysical dependences into the quantities we want to compute. We will discuss 
this in more detail in section IV. 

Using ( |3.22| ) and ( |3.23| ) in ( p. 211 ) we write the 0{as) 2 to 3 body cross section as follows: 



[q{pi)q{p2) ^ W^g] 



^qq {finite) + ^qq (col+) + ^qq {col-) + ^qq {soft) + ^l^j 



(3.24) 



with 



a. 



qq {finite) 



1 1 

»i 



2-i0---4„-i 



X 



X 



1 



dy 



TT S 



1-2/ 



P{s) 



dx P{sx) 



X 



xo 



yo 



1+2/ 



yo. 



dcos9[ d9'^ F'''^{s,x,y, cos e[, 9'^) 
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1 1 



qq (coZ±) 



AN^ 2s 

rl 



(47r 



>e-2 



X 



X 



l-xJxo I X J 

COS 9[ Sin'^e ^/^ ^gg- M±) ^/ 



a. 



qq (soft) 



1 1 

4iVf 27 
"1 



(4vr; 



.-2 r(i - e) s 



-l-e 



TT 



(l-Xo)-'^$(s) 



r(l - 2e) 27r 2e^ 



X 



(3.25) 



where 



F'?'?(™'±)(s,x,cose;) = = ±l,cos^i,cos^^) 

F««"(^°^*)(s,cos^i) = = l,i/,cos^'i,cos^^) . (3.26) 

To compute the quantities in ( ^.26| ) we first note that the following relations hold for y = 1: 



Pi - k 



xpi 

(1 - x)^i 



(3.27) 



Using these relations in (|3.16| )- (|3.18| ) and noting that in the limit y 1 (—1) only Mj^ 



) contributes to F'^'^ in (|3.2(]| ) we obtain the covariant expression: 

1 + a;^ — e(l — x^'^'^ 



F''«"(^°'+)(s,x,cose;) = SsglCf 



X 



(3.28) 



We implicitly include in 
computation yields: 



F''^"(™'-)(s,a;,cos^;) = 8s^|a 



I) an overall factor of C((5i,<52,(l — An analogous 

{I + x"" - e{l - xf) 



X 



Y.\M^^^^\xs,h-) 



(3.29) 



In the latter we implicitly include an overall factor of C{Qi,Q2, (1 — x)p2)- The collinear 
limits of F'^^ have thus been reduced to 2 to 2 body squared matrix elements. We checked 
that the same expressions for the collinear limits of F'^^ are obtained when the sum of gluon 
polarizations is taken in the axial gauge. 

Now we obtain the soft residues of each of the contributing terms in the squared matrix 
element: 



\imA{p,.k){p2-k)J2\Mf, 

X— >1 ^ — ' I 

lim4(pi-fc)(p2-A:)Ekf 
lim 4(pi ■ k){p2 ■ k) MZM*f + c.c. 



AsglCp{l + yf |m^^"^^^(3, ¥"f') 
WsCF{\~y^)Y.W'~ 



(3.30) 
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After summing the above three contributions we obtain: 

piq (soft) ^^g^/) _ WsglCp E Im'^'^^^s, b'°f'] 



(3.31) 



This contribution contains an imphcit factor of C(Qi, Q2, 0). Note that the dependence on 
y cancels after all terms in ( ^.301 ) are added together. We have checked that the same result 
is obtained if the gluon polarization sum is taken in the axial gauge. In ( p.28|) , (|3.29|) and 
(|3.31| ) we define: 



b- 



b{s, x,y = 1, cos 9[) = (1 + cos6'^) 



sx 



b{s, x,y = —1, cos 9[) = — (1 + cos6''5^^ 
b{s,x= l,y,cose[) = (1 + cos6'i) . 



(3.32) 



Note that in the soft limit the variable cos d[ is equivalent to cos of a 2 to 2 body kinematics. 

Noting that the squared matrix elements in ( |3.28D -( p.31D are of the 2 to 2 body type we 
can now rewrite the soft and collinear terms in a more convenient way: 



p 

"qq [soft) 



^C, -2V - 
7r V e 



+ 21n(l-xo) 



+ |ln4 + 21n2(l-a;o) 



+2 ln(l - xo) In ^ + 2C(2) - 4 ) (t(°) [q{pi)q{p: 



a, 



as 



qq (coi+) 



27re 



I dx\cJl + x')\-^] +ecJln(^J^]( 

e J Pis) [ ^ ^ ll-x}xo \ \syoJ ^ 

ln(l — x) 



1 X J I'o 



- 2 1 + 



1 — X 



\ + X' 

+ x-l\ 



a 



qq {col-) 



dx\cF{l+x^)\^\ +eC^flnf^^ 
eJp{s) [ ^ ^ U-xJa.o \ \syoJ 

- 2 (1 + x') 



{ln(l — X 
1 -X 



1 + x^ 
+ x- 1 



xo 



xcr^°^ [g(pi)g(xp2) ^ Wj] 



In the previous formulae we neglected terms of 0(e) and we used 



(3.33) 



1 



1 



2e2 2e 



12- 



4 /i2 



(3.34) 



We also made the replacement = AnasfJ^'^'' ■ 

The 0{as) corrections to the 2 to 2 body partonic cross section for qq W'y were 
computed in [0 and they are given by: 

[q{pi)q{p2) ^ Wj] = ^C^a^'^ [q{pi)q{p2) ^ W^] 2V 
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where 



^1 (2 (s-M^-b)- b) (2Fi(g, 5) - F2{s, b)) 

xC(gi,g2,o) 

(3.35) 



F,{s,b) = F{-b,-{s-M^-b),s,M^) 
F2(s,6) = Fi-is-M^-b),-b,s,M^) 

Cw ^ (3-36) 



and Gf is the Fermi couphng constant, b was defined in (|3.6| ) and F{ti,t2, s, M^) was defined 
in O of i. 

Using (pD, dpI), dpI) and (lOSp in the RHS of (|18|) we can write the 0{as) hard 
scattering cross section on the LHS of (p.l8|) after performing an integration over Qi and 

a^^^ [q{pi)q{p2) W-f] + a^^^ [q[pi)q{j>2) W^g] = 

^£ (finite) + (^qq (soft) + ^^^^^^ [?(Pl)^(P2) ^ W -f] 

p , 



+<a + ^td^ PqqiM'^ [q{Pi)q{vp2) -> W^] 



qq (col^) -r 

P P 
^qq (finite) + ^qq (SV) 



27r 7p(s) 



/' dx (l + x^)ln(^)(-^j -2(1 + 



/II I - V- ■ ^ ^^T^ — ^ I* +a;- 1 
s?/o / U - a; J ^0 ^ M 1 - a; 



ln(l — a;) 



(3.37) 



where we have defined the soft-plus-virtual contributions: 



p _ p 

^qq (SV) - "qq (soft) 



^Cf ( ^ In 4 + 2 ln2(l - xq) + 2 ln(l - Xo) In 4 

TT \2 /i^ /i^ 



xaW [g(pi)g(p2) - 1^7] + V^'^^^ 
(2 (s - - 6) - 6) (2Fi(s, 6) - ^^(s, 6)) 



J-i s - 

(3.38) 

Summarizing all the contributions in the incoming qq partonic channel, we have 



J DQ.DQ, [ dq.DQ^ J ^^^^'^ ^^^'^ " ^^^'^ ^^^'^ 
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+ dTf,,{r) j DQ.Dq, 



(3.39) 



with pi = TiPi, p2 = T2P2 and q2 = Q2/T. Since there are no singularities left all the 
necessary squared matrix elements in (|3.39|) can now be safely evaluated in n = 4 dimensions. 
On the RHS of ( |3.39D we still need the following squared matrix elements: 



^ |m«^"^^^(s, b) 



2(0) 



—irNcC^a 
y 



(2(^ 



b)-b 



X 



b{s 



M^-b) + Us 



b{s 



2(1) 



'b^ + {s-M^-by + 2sM^ 



w 



b{s 



M^-b) 



(3.40) 



where b = b/r. The integrand of the last term in ( p.39| ) carries an implicit factor of 
C{Qi,Tq2^{l — T)q2)- In (|3.38|) and (|3.40|) the invariant b is evaluated in the unprimed 
frame as defined in (f3.6D. 

2 



The expression for Yj 



in n = 4 needed when evaluating a 



qq {finite) 



is too 



long to be presented here, but it may be obtained upon request. We note here that 75 is 
never needed in n 7^ 4 dimensions: we obtained the cancellation of singularities before fully 
computing any squared matrix element where we had to explicitly evaluate 75 and whatever 
remained after this cancellation could be safely computed in n = 4 dimensions. 

Since we have integral expressions of all quantities on the RHS of p.39| ) in terms of the 
variables x, 9'^ and ^2 which define all the independent invariants of the system we can 
compute these integrals using numerical Montecarlo techniques and histogram any physical 
variable of interest that can be expressed in terms of these invariants. We will reexamine 
these issues in more detail in sections IV and V. 



III.C The qg CHANNEL. 

In this channel the singularities are not both of initial state (type /), as in the case of qq 
channel, but we have now one piece in the initial state (type J) and another in the final state 
(type //), as shown in Fig.2. In this case it is thus more convenient to integrate each term 
separately and write for the total partonic cross section in this channel: 

Wi)g{P2) W-iq] = a^'^ [q{pi)g{p2) W-fq] + a^'^^ [q{pi)g{p2) ^ W-fq] 

+ [<l{Pi)9{P2) W-iq] 

(3.41) 
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where the first two terms on the RHS of ( p.41|) contain only the partial squared matrix el- 
ements J2 l^ibf {Mfff respectively, while the third term contains J2 l^nif ^ind all 
interference terms of the squared matrix element. If the gluon and the photon are summed 
over physical polarizations only then the third term on the RHS of ( |3.41| ) is free of singular- 
ities while the first two terms contain only collinear singularities. 

If we define the same kinematics in c^^^^^g as we did with the qq channel, the collinear 
pieces of integrand contain a factor [1 + y)~^. To isolate the singularity in this term it is 
therefore enough to define the non-singular function: 

Ff'is,x,y,cose[,e',) ^ -2p2-k^\Mr,\' (3.42) 

so that the first term on the RHS of ( p.41| ) may be written: 

o^'^"^ [(l{Pi)9{P2) W-^q] = 



(47r)'"^ — ^ — / dx ^{sx){l - X 



SN^Cpil - e) 2s ' ' r(l - 2e) n Jp{s) 

C dy{l-y)-'{l + y)-^~' dcos9[ sin-^' 9[ F dO'^ sin'^^^^ 
7-1 J-i Jo 



We now rewrite the factor [1 +y) ^ as a distribution: 



(3.43) 



[l + y)-'-' - -^6{l + y) +\—-\ +0(6). (3.44) 



yo 



Using ( I3.44D in ( p.43| ) we can write the corresponding contribution to the 0{as) 2 to 3 body 
partonic cross section as follows: 



cr""'"^'^ [q{pi)9{p2) - W^q] = a^^^^, + + 0(e) , (3.45) 



where 



X r de'^ Ff (s,x,y,cos^'i,^^) 



x{l-x)-^' j'^dcose[ sm-^'e[ Ff ^'=°'-)(s,x,cos^;) 



(3.46) 



and 



(s, X, cos e[) = Ff (s, x,y = -1, cos e[, 9'^) . (3.47) 
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To compute the limit on the RHS of ( p.47|) we write, as we did in (|3.16|) and ( p.l7|) : 



9l 



(2p2 ■ k) 



xM* [q (pi) q (p2 - k, a') ^ WiQMQ2)] 



{3A8) 



where now in the qg center of mass frame: 



T>aa 



-4p2 ■ k 



7o 



"2-n ^ {l-x){l-yy 



-(l-a;)(l + i/)^i| 



2 



((l-x)(l-y)-l)^ 



(3.49) 



From the above equation it is now clear that there will be no soft quark singularity coming 
from the type la squared matrix element. Using ( p.49|) , ( p.47|) and (|3.42|) we obtain: 



-,qg {col-) 



s,x,cos9[) = 2g^CF 



(2x(l - x) - 1 + e) 



X 



(3.50) 



The latter expression contains an implicit factor of C((5i,<52, (1 — x)p2) to account for ex- 
perimental cuts. Using ( p.50|) in (p.46|) we can rewrite the coUinear contribution (neglecting 
terms of 0(e)): 



a. 



qg,col- 



as 



27re Jp{s) 



dx 



1 



x{l — x) — - + e < In ( x{l — x) — - ) + ln(l — x) 



1 



+x(l - x) [1 - 21n(l -x)] 



[q{pi)q{xp2) W^r] ■ 



(3.51) 



To treat the type // term on the RHS of (|3.41|) it is convenient to rotate the reference 
frame so that in the Wj center of mass frame we have: 



Pi 
k' 

Q2 



p[Q{l,0,...,0,smip',cosip') 
P2,o (I'O' ... ,0, - sin x',cos x' 
A;;(1,0,... ,0,0,1) 



Q'l 



sin (f)[ sin (p'^ cos (f)'^, — sin cos — cos (p'^ 



Q[ (1, . . . , sin sin 02 cos 03, sin (f)[ cos 02, cos 0^) . 



(3.52) 



We can thus write for the 2 to 3 body integral: 



DQIDQ2 



DQIDQ2 



[q{pi)9{p2) ^ W{QMQ2)q] 
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1 



(47r) 



^2 r(l-e) sl-%2-2e. 



8N^Cf{1 - e) 2s ' ' r(l - 2e) 27r 



/ dz ^ {s[l + I3{s){z - 1)]) 



X 



x^|M|f(s,a,6,c,c?)|' 



sm 



-2e i/ 



(3.53) 



with 



cosx 



1 — X — y{l + x) 
1 + X — y{l — x) 



V = cos I 



1 + 



X — 1 



(3.54) 



The exphcit form of the invariants in the squared matrix element in terms of the new inte- 
gration angles is given as follows: 

b = 2pi-Q2 = 2p[q Q[ (1 — sin-i/^'sin^'^ COS02 — cosV^'cos^'i) 
c = 2A;-Q2 = 2A;o|g;|(l-cos0;) 

d = 2p2 ■ Q2 = 2^2 Q'l (1 + sin x' sin 4>'i cos 02 — cos x' cos . 

(3.55) 

The rest of the invariants and variables remain as defined in ( |3.11D -( p.l3| ). To isolate the 
singularity in the type II squared matrix element it is enough to define: 



F?f{s,z,y,v,<P'2) = -2k . Q^Y^lMfff 
so that now the second term on the RHS of ( |3.41D may be written: 

[q{pi)9{p2) - W^q] = 

1 1 (4^^-2 r(l - 6) 2^^ 1 



(3.56) 



SN^Cpil -e)2s' ' r(l - 2e) tt r(l - e) 167r 
f'^ dz\z{\ - z)\-^^ [1 + (z - V)(i{s)X W dy (1 - y') ' 



X 



X / dv{l + v)-'{l-v) 



-l-e 



2 sin '02 F^f {s,z,y,v,(j)2) 



We rewrite the factor (1 — t>) as a distribution: 



(1-^) 



^0 



■5(1-^;) + 



1-v 



0(e) 



with < fo < 2, so that ( ^.571 ) may be rewritten as follows: 

C^PJIW [q{p-^)g{p2) W-fq] = (T^fjinUe + ^qg,Ll+ + ^(e) 



(3.57) 



(3.58) 



(3.59) 
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where 



a. 



-P,ii 

qg, finite 



1 1 



(Any' Pis) £dz f^dy j\dv[-^ 



X 



a. 



qg,col+ 



%N^CF{l-e)2s^ ' r(l-e) vr e\v^ 



^{s)fi-\s 



X dz\z{\ - z)\-^^ [1 + - \)fi{s)X J ^ dy (l - y') 



qg {col+) 



II 



s,z,y) 



and 



Fff^-^^\s,z,y) = F^f{s,z,y,v = lA'2) ■ 
The type // squared matrix element is given by: 



(3.60) 
(3.61) 



E 



{2k ■ Q2 



^R^g'ii E M [q{pi)9{p2) ^ W{Qi)q{Q2 + A;, «)] 
xM* [g(pi)(?(p2) - H^(gi)g(Q2 + k,a')] 



(3.62) 



with Cg the charge of the outgoing quark. In the qg center of mass frame we have: 



RqgJI = -4fc-Q2 



7o 



'2~n (l-z)(l + v) \ / il-z){l + v) 



2 2z 

' {i-z){i-vy 

2z 



(3.63) 

where we have summed over physical polarizations of the outgoing photon in the covariant 
gauge. Again we note that there will be no singularities in the soft quark limit, that is when 
z —>■ 1. For the coUinear limit of Fjf we obtain: 



M^O^Wql^^JL 



with 



hji = b{s,z,y,v = 1) = -p{s)z{l -y) . 



(3.64) 



(3.65) 



Remember that ( p. 641 ) has an implicit factor of C{Qi, Q2, {1 — z)Q2/ z). Using ( p. 641 ) in ( |3.60| ) 
we can rewrite the coUinear contribution (neglecting terms of 0(e):) 



a, 



qg,col+ 



a 



27re Jo 



dz 



„2l + (l-^)' .2 



1 + 11-^ A /2^- 

z \ \svo ^ 
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+ ln 



l + {z-l)l3{s 



^2(1-^)2/32(5 



q{pi)9{P2) ^ ( — 



(3.66) 



In ( p.66|) we made the replacement = Anafi'^'^. 

Using (lOTD , (|3^),(g), (|33ll) , (CT),(p:60D and ( M ) in (F^) we obtain for the 
0(0:5) hard scattering cross section: 



-gfjinite + ^ J^^^^ rfx |- + (x^ + (1 - x)' 

Xd"^"^ [g(pi)g(xp2) W'j] 



a. 



ln(l - x) - - ( 1 + In 



'2/i2 



P,ii 



a 



'^^qg, finite + 27r^'' 



dz < Z 



l + (l-z) 



In 



+ In 



z'^il- zY^-^^s) 



X / DQ,Dq2——— [q{pi)9{p2) ^ W{Q^)q{q2)] 
J DQiDq2 



(3.67) 



with q2 = Q2/T. We have thus cancelled all singularities and we can now summarize for all 
the contributions in the incoming qg partonic channel: 



[p(Pi) p{P2) ^ W{Qr) 7(Q2) X] 



1 /•! 

dri / dT2 {fgp{Ti)fgp{T2) [q{pi)9{p2) ^ W-fq] 
Jo ^ 

+ dTf,,{T) J DQ,Dq2^^-^^[q{p,)9{P2) ^W{QMQ2)]] 

+ {p ^ P,Ti ^ T2,Pl ^ P2) } 

(3.68) 

with pi = TiPi, p2 = T2P2- In (|3.67|) and ( p.68|) we need the following squared matrix 
element: 



V s(s -M^-b) 



(3.69) 



where b = b/r. (|3.69| ) carries an implicit factor of C{Qi, rq2, (1 — r)g2)- The pieces of squared 
matrix element needed in the last term in ( p.67|) are too long to be presented here but they 
may be obtained upon request. 
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III.D The gq channel. 

The treatment of this channel is analogous to the qg channel. We can again decompose the 
partonic cross section as follows: 

[9{Pi)q{P2) W-fq] = a^'^ [g{pi)q{p2) W-fq] + a^'^^ [g{pi)q{p2) W-fq] 

+ 5-^ [9{Pi)q{P2) W-fq] 

(3.70) 

where the non-singular term [g{pi)q{p2) W'-fq] contains all interference pieces of the 

squared matrix element and also Mfjj . The other terms in ( p.70|) are decomposed as 
follows: 



(^^'^^^^ [9{Pl)q{P2) W-fq] = (Tg^jinite + ^^icol+ + ^(e) 

^P,nm [g{p^)q{p^) W-iq] = (Tg^jinUe + ^giLi+ + ^(e) 



(3.71) 



with 



P,i 

^gq, finite 



^gq,col+ 



(477) 



Pis) 



dx P{sx) / dy 



d cos 9[ 



yo 



X / de', Ff {s,x,y, cose[,e',] 



as 



27re Jp{s) 



dx 



x{l — a;) — - + e < In 



x(l — x) — 



Infl — x] 



+x{l-x) [l-21n(l-x)] 



cr 



(0) 



[qixpi)q{p2) -> W'j] 



a 



-PJi 

gqjinite 



1 1 



(47r)~^/5(s) dz J^^dy f^dv^^j^ 



V 



VQ 



X 



^2 Fff (s, z, y, f , 02) 



p,ii 

^gq,col+ 



a 



2TTe Jo 



dz 



.i + n-z) 



1 + (1-^) 



In 



+ ln 



' i + {z-i)f3{sy 
^z^{i-zyp^{s) 



X / DQ^DQ,-^^^^ \9{pi)q{P2) - W^(Qi)g (— 
J DQ1DQ2 I \ z 



(3.72) 



In ( |3.72D we used the following non-singular functions: 

Ff (s,x,y,cos^;,e^) ^ -2pi-A;^|Mf 

F//(s,z,i/,t;,0^) ^ -2fc-g2E|^/ir ^3-^3) 
Using ( |3.7(J| )-( pr7^ ) in ( ^.3(J| ) we obtain the cancellation of all singularities in this channel: 

(^^^^ [^(^1)9(^2) W-iq] = 
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lp{s) 

xcr^°^ [q{xpi)q{p2) —> W^] 

'1 + (1-^)2 



ln(l -x) - - ( 1 + ln 



a 



' gq, finite 27]- « 



Inf^l +ln 



X / Dg^Dg^— — — [9{pi)q{P2) ^ W^(Qi)g(g2)] 

J DQiDq2 



[9iPi)qiP2) W-fq] . 
We can now summarize all contributions in the incoming gq partonic channel: 



(3.74) 



/ DQiDQ2j2 



b(Pi) p{P2) ^ W{Q,) 7(g2) X] 



DQ^DQ2, 

dn I dT2 \fgpin)fgp{T2)ia^^^ [g {pi ) q{p2 ) ^ 1^7?] 



1 rl 



JO 



+ J^dTf,,{r) J DQ,Dq2^^^ [g{pi)q{p2) ^ W^(Qi)g(g2)] 



(3.75) 

with pi = TiPi, p2 = T2P2 and q2 = Q2/t- In (t3.74|) and (|3.75|) we still need the squared 
matrix element: 



2(1) .fi,.^^. fs^ + P-2{s-M^-b)M? 



w 



sb 



(3.76) 



with b = b/r. The comments after ( |3.69| ) apply here too. 
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IV THE THREE SCENARIOS AND THE EXPERI- 
MENTAL CUTS. 

IV.A 2 BODY INCLUSIVE PRODUCTION OF AND 7. 

In this scenario ("2 body inclusive scenario") one does not tag the outgoing jet, so it will 
include events with zero and 1 outgoing jet. We may define this scenario requiring the 
following conditions for the outgoing particles 



, cos 6^,^ 


< 


cos( 


Pt-y, Ptyy 


> 


Pt- 


Rw,-y 


> 


R- 



{Rjet,-y < R ) =^ (■S(jet,7) < S ) 
{Rjet,W < R') =^ iSijet,W) < ■5") (4-1) 

where we call 6i the angle between the incoming proton axis and the axis of the outgoing 
particle i; Pti is the transverse momentum of particle i. Rij is the cone size between a pair of 
outgoing particles: Rij = {Aijri*y + (Aij^)^ with the pseudorapidity rj* = (1/2) ln[(l + 
cos^^)/(l — cos^^)] and (f) the azimuthal angle; S(^jet,w) = Ejet/Ew is the "shadowing ratio" 
between the untagged jet and the W boson. The last two conditions in ( [4.1| ) discard events 
where the jet being too close to the W or the photon is at the same time of comparable 
energy so that it would "shadow" one of the two tagged particles making it undetectable. 
For this purpose we check the cone size Rjet,'y {Rjet,w) and if this is less than R~ we keep 
the event only when S(jet,7) (s(jet,VK)) is less than s~, setting the differential cross section to 
zero otherwise. The quantities 6~,Pt~,R~,s~ are constants related to the acceptance and 
resolution of the detector. All the quantities are defined in the proton-antiproton center of 
mass frame. 



IV.B Production of W+ and 7 with 1 jet. 

Here one detects three outgoing particles, namely W^, 7 and 1 jet. We call this the "1 jet 
scenario" and we define it by imposing the following conditions: 



cos 9w\, cos 6*^6* 


< 


cos(i 


Pt^, Pt\Y, Ptjet 


> 


pr 


Rw,'y 


> 


R- 


Rjetj-y 


> 


R- 


Rjet,W 


> 


R- . 



(4.2) 



IV. C Production of and 7 with jets. 



In this scenario ("0 jet scenario") we select events where the and 7 are detected but no 
outgoing jet is detected. This includes 2 to 2 body events and 2 to 3 body events where the 
outgoing jet has a small angle with respect to the beam, a small transverse momentum or it 
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is "shadowed" by the photon or the W so that it remains undetected. We may define this 
scenario requiring the following conditions for the outgoing particles 

I cos 9^ 



, COS 9w 


< 


cos{9 


Pt^, Ptw 


> 


pr 


Rwrf 


> 


R- 



iRjet,-y < R ) =^ (■5(iet,7) < S ) 
{Rjet,W < R ) =^ {sUet,W) < ^ ) 

(I cos^.etl > cos(^-)) or {Pt^.t <Pr). (4.3) 

IV.D General remarks. 

We note that the second and third scenarios are complementary, in the sense that an event in 
the first scenario falls in either of the last two. In other words, we may obtain the histograms 
of the jet production scenario by subtracting the histograms of the 1 jet scenario from the 
corresponding histograms for the 2 body inclusive scenario. 

To implement the three experimental cut functions C{Qi, Q2, k) which define each of the 
scenarios in A,B and C all quantities involved in the above conditions have to be defined 
in terms of the partonic invariants that are used in the integrands of the corresponding cross 
section formulae. In the proton-antiproton center of mass frame we have: 



cos 9^ 



Ew 
cos 9w 



Ejet 



Pl-Q2 + P2- Q2 

Pl-Q2~P2- Q2 
Pl-Q2 + P2- Q2 
Pl-Ql+P2- Ql 



Pi-Qi-P2- Ql 



'{Pi ■Qi + P2- Ql? - SM^ 
Pi-k + P2-k 



P,-k-P2-k 

cos9iet = ; r- 4.4 

^ Pi-k + P2-k ^ ' 

Pi and P2 represent the proton and antiproton momenta respectively; they must be appro- 
priately expressed in terms of the incoming parton momenta pi,p2 and their momentum 
fractions ti,T2 in all the cross section formulae, = 2Pi ■ P2 is the proton-antiproton 
center of mass energy. Qi, Q2 and k are the momenta of the W boson, the photon and the 
jet respectively. The rest of the quantities needed can be computed using the ones in (|4.4|). 

When we replaced the divergent factors (1 — (l±?/)~^~^ and {l—v)~^~'^ in section 
III with distributions the resulting equations remained valid as long as the variables x, y and 
V were integrated over their whole range. The energy of the outgoing jet in the incoming 
parton-parton center of mass frame is linearly related to the variable x (see ( |3.10| ),) so it 
is in principle not a physical quantity unless x < Xq, in which case the symbol ~ can be 
replaced by = in the corresponding distribution in (|3.22| ). Similarly, the angle between the 
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outgoing jet and the beam in the parton-parton frame is related to y and the angle between 
the outgoing jet and the photon is related to f , so these quantities are not physical either, 
unless the variables y and v fall inside the ranges where we can replace ~ for = in the 
corresponding distributions. 

According to the above observations we shoudn't have any trouble in the 2 body inclusive 
and in the jet scenarios, since in these cases the outgoing jet is not being tagged so the 
unphysical variables are not "observed", but they are rather integrated over their whole 
range. However, in the 1 jet scenario the energy and angles of the jet are observed and these 
are directly related to the variables x, y and v. According to the way we defined the 1 jet 
scenario in ( [4.2| ) the outgoing jet is never allowed to be soft or coUinear to the beams or the 
outgoing photon so the subtraction of divergences will never be active. With this in mind 
we can easilychoose the parameters xq, yo and vq in our Monte Carlo in such a way that the 
sampled ranges of x, y and v always fall inside the regions where ~ may be replaced for = 
in ( p.22| ), ( |3.44|) and ( p.58|) . To accomplish this we can just take Xq = 1 and yo = Vq = 0. 



The experimental cut function C((5i, k) will a 

ccordingly set to zero all the terms containing ill defined logarithms. 
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V The numerical implementation. 

When numerically implementing the "generalized plus" distributions defined in ( |3.23D to 
compute total cross sections the second terms on the RHS of ( |3.23| ) are finite when the 
soft or collinear limits are approached. However, when we produce histograms of single or 
double differential cross sections it is necessary to split the second terms on the RHS of these 
definitions into two parts, as we will explain next. For the case of the x integration we have: 



' XQ 1 — X Jxo 1 — X Jxo 1 — X 

The first term on the RHS of ( |5.1|) is naturally histogrammed using 2 to 3 body kinemat- 



ics. The soft pieces that resulted from the expansion in ( |3.22| ) were added to other 2 to 
2 body contributions in order to cancel singularities so the remaining pieces are naturally 
histogrammed using 2 to 2 body kinematics. This means that in order to keep consistency 
in our computation we have to histogram the second term on the RHS of ( |5.1|) -which is the 
term that compensates for the soft singular terms in (p.22|) - using 2 to 2 body kinematics as 



well. It is thus clear that a consistent histogramming cannot be achieved in a simple way 
without splitting the LHS of (|5.1| ). In doing so we introduce logarithmic singularities in each 
of the terms on the RHS of that cancel each other only after summing both contribu- 
tions bin by bin. To control the numerical cancellations we introduce small adimensional 
cuts Ax, Ay and A^, in the lower or upper limits of the corresponding integrals. A first order 
estimate of the error introduced by the cuts along with the requirement of good numerical 
convergence will help us find the best values for these parameters. 

In what follows we will rewrite the partonic hard scattering cross sections for each chan- 
nel taking into account the A parameters introduced above. The contribution of each of 
these terms to the hadronic cross section is obtained after multiplying by the corresponding 
experimental cut function, convoluting with parton densities (see section II) and adding the 
corresponding "inverted channels" (i.e. the ones obtained by interchange of the incoming 
partons.) Numerical results for each of these hadronic contributions are presented in the 
following paper 

For the qq hard scattering channel cross section needed in ( p.39| ) we have: 



Born I P I I I I P I ^ I ^ 

^gg ~ ^qq ^ ^gg [SV) + ^la + ^Ih + Cr/,4 + O^qq {finite) + '^gg {Brems) + CTqq {error) 



(5.2) 



where 



a 



.p 



O'la — 0"/a,l + 0"/a,2 + 0"/a,3 



qq {finite) 



0"/,l,l,a + ^/,l,2,a + '^f,l,3,a + f,l,l,b + ^/,1,2,6 + ^/,l,3,f) 
+<7/,2,l,a + Cr/,2,2,a + <^/,2,3,a + f,2,l,b + Cr/,2,2,6 + f,2,3,b + 0"/,3 



^qq {error) ^la,error ~l~ ^ Ib,error ~l~ f,l,error,a ~l~ ^f,l,error,b ~l~ '^f,2,error,a 

f,2,error,b ~l~ '^f, error 

(5.3) 



34 



and 



Born 



CTlaA 



Cla,3 = 



error — 



crib. 



error — 



Jo ^ I 



2(0) 



Cqq,l 



f 

Jp(s 



I3{sx) 



J p{s) X J —I I 



2(0) 



2 yp(s) x 1 1 - x 



ln(^^^ -21n(l-x) 



X J^^dcose[ Yl \M'^^^^^ixs, xb+) 



2(0) 



xo 



X 



1 + x^ 
1-x 



+ 2 - 2 In 1 



ln( ^1 -21n(l-x) 

2(0) 



2 I Uyo 



X 



d_ 

dx 



p{sx) 



X 



(l + x^) J'^ d COS 9[ J2 \M'''^^^{xs, xb+) 



2(0) 



x=l 



+0(A^lnA, 

Cqq,l 



2 J p(s) X J —I ' 

hi 



2(0) 



X 



2 Jp(s) X \1 — X J 
J^^d COS 9[ 6") 

2 Jxo X \1 — X J 



'2/^2^ 



21n(l -x) 



2(0) 



Inf ^ ) -21n(l-x) 

2(0) 



'2^2 



'"'"'A^ (In I ^ 1 + 2-21nA, 



X 



d_ 
dx 



P{sx) 



X 



(1 + x^) y d cos Y |m««"^'^^(xs, r 



2(0) 



+0(A^lnA.) 

(77,4 = 2Cqq- 1 /5(s) In 



1-Xo 



In f - In [(1 - xo)A, 
\syoJ 



X l^^d cos 9[Y\M''^^^is^b'°f''^ 



2(0) 



^0 , /3(sx) /-l-'W) 

ax / a|/ 

p(s) 1 — a; J-i 



dcos^; r d02 F'i%s,x,y,cos9[,02) 

1 — V J -I JQ 
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(^f,l,error,a 



<^/,2,3,6 



^f,2,error,b 



•^0 /5(sx) 



1 /•! 



1 - a; 7i_j^o 

.2 



1 - 2/ J-i 
1 + x' 



4 V^y/ J pis) I — X \ X 
s 7p(s) 1 - X y-1 ^0 \ oy ^ 

I a2\ 



dcose[ r de' F^^s, x,y, COS e'e') 

Jo 

i \ 



^gg,2 r 

s Jp{s) ^'^ 1- X J-i+yo l + y 

Jp{s) 1 — X J-i+Ay 1 + y J-i Jo 



y=l 



' V yJ 

r° dx^^ [' dy^— C dcos9[ F dO' F'i'^(s,x,y,i 

S Jp{s) 1 — X J-l+yo I + y J-1 Jo 

' - ^ ™- e[A) 



A A) 



4 \AyJ Jp{s) 1-x 



'1 + x^ 

X 

2(0) 



(^f,l,error,b — 

^/,2,l,a = 

C^/,2,2,a = 

^/,2,3,a = 



J^^d cos e[ \mi'^^^^{xs, h-)^ 
^A, fda;^ tdcose'^Tde', {■ 

S Jp{s) 1 — X J-1 Jo \ 

+0(Al) 



f dF'i%s,x,y, cose[,e'^) 



+o{Al) 



9y 



■j=-i 



/' dcose[ r do'^ F'^%s,x,y, COS e^e'^) 

s Jxo 1 — x J-i i- — y J-i Jo 

/-i-A. ^^^(sx) f^-''- ^yj_ ^^^g^, r^^, F'i%s,x,y,cos0[,( 

S Jxo 1 — X Jl-yo ^ — y "'-1 Jo 

\\/ Jxo 1 - X \ X J 

dcosO' V lM««^^Tr.T.s..Tft+l 



error.a — 



4 \I\yJ Jxo i. — X \ X J 



a:,y,cosg;,g^) 



9y 



+0(AJ) 

■ J^^ d cos e[ £ de'2 Fi'i{s, X, y, cos e[, 9'^) 
- C dcos9\ r dO'n F'iHs.x.v.cose[,e'^) 



^r^dx^-^r dy^[\ 

S Jxo 1 — X J-l+yo 1 + y J-1 

^J_{SX)_ pl+.o ^^^^ ^1 

Jxc 



1 — X J-l+Ay 1 + y "'-I 
ixo 1 — X 



S Jxo 1 — '• 



■1+2/0 1 /""■ 

dy-^— I dcose[J de'2 Fi%s,x,y, 



X 



a; 

2(0) 



s ^ Jxo 1-X J-1 Vo V 9y 



u=-l 
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-2C,,- 1 (3{s) In 



1-3^0 

A. 



ln(— jy c/cos^;^|M«^"^^^(s,6^"^*; 



2(0) 



+ 



1 



1 + y 



yo. 



dcos6[ / 



dx 



{sx)F'^%s,x,y,cos9[,9',)] 



x=l 



(Tgg (Brems) = drf^g^r) J ^Q^^^'^ JJQ^JJ^^ [(l{Pl)<l{P2) ^ W{Qi)g{q2)] 



2(1) 



(5.4) 



The invariant h in the unprimed frame was defined in ( |3.6D . In ( |5.4| ) appropriate experimental 
cut functions are imphcit in each of the corresponding integrands. We have introduced the 
constants: 



C. 



C, 



qq,3 



1 1 1 

4iV|2ll67r2 
111 

4iVf 2^210^ 
111 

4iVf 2^16^ ■ 



asCi. 



(5.5) 



Now we rewrite the hard scattering cross section for the qg channel needed in ( p. 681) : 



cr. 



PJ 



'^qgjinite ~^ '^qg ~^ '^qgjinite ~^ '~^qg' ^qg '^qg{Brems) + (^qg{error) 



PJI 



ILcol I ~P 



qg 



(5.6) 



where 



P,i 

^qg, finite 

P,II 
'-^qg, finite 

^qg(error) 



(^IgJA + ^lgj.2 + ^lg,f,Z 

^qgj,l + ^qgj,2 + CTggJ.S 

I I // 

qg, error ' ^ qg, error 



(5.7) 



and 



a 



19 J ^ 
.1 

q9,f,2 
I 



Cq„3 / dx (3{sx) / dy / dcosO'^ / d62 Ff^ {s,x,y,cos9[,92] 

' Jp{s) J-l+yo l + y J~l Jo 



^qg^ 



Pis) 

1 

Pis) 



dx P{sx) 



dy- 



dcos9[ d92 Fj^ {s,x,y, cos 9[, 62] 
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a, 



qg,error 



= Coo 3 A„ / ax pisx) / d costal / dOn 

J Pis) J-1 Jo \ ay 



y=~l 



a. 



^I,col 
19 



J pis) X [2 ^ 



a, 



a, 



a 



II 

QgJA 
II 

II 

19,1,3 



-Cqg,3 f^i^) dz j^^dy 



rl n 
-Cgg,3 P{s) I dz dy 



1-iln, , 

2 2 \syoJ 

1-DO 1 

dv 

-1 1 — f JO 

1-A„ 1 



ln(l — x) 
1^ d cos 9[ ^ |m''^"^^^(xs, b- 



2(0) 



dv- 



1 Jl-vo 1 — V Jo 



^2 Ffi is,z,y,v,(f)2) 
02 Fii (-5, z, y, V, 02) 



2(1) 



qg, error 



„II,col 
"^19 



c,,,KPis) tdzrdyrd^'J'^-^'^'^y^'^^'^^ 



1 Jo 



Cqg,2 el Pis) 

X dzlz — 



'l + {l-z) 



In I 

2(1) 



SVo^ 



dv 



+ In 



+ 0(AS) 



v=l 



'i + iz-msY 

^z^{i-zyp^{s)^ 



^19 



(TqgiBrems) = drf^gir) J ^Q^^^^ jjQ^jj^^ ^(^1)^(^2) W{Qi)q{q2)] 



Cqg,i P{s) drf^q{r) dcosOiY, 



Mi9^^i{s, -) 
r 



2(1) 



(5. 



The comments after ( |5.4| ) are also valid here. In ( |5.8D we have introduced the constants: 

1 1 1 



C'g3,l 
^19,3 



as 



a 



8N^ 2s 167r2 
1 1 1 

8N^Cf2^32^ 
1 1 1 

8N^Cf2~s¥^ 
1 11 

8NfC^2^1Q^' 



a 



(5.9) 



Finally, the hard scattering cross section in the gq channel needed in ( |3.75|) may be 
rewritten: 



J,col , P,ii 



II,col , ~p 



^91 ~ ^gqjinite + ^g'q '^gq, finite "I" '^gq' "I" '^gq "I" ^gqiBrems) + (^gqierror) 



(5.10) 
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P,i 

^gq, finite 

P,II 
^ gq, finite 

^ gq{error) 



<^gq,f,l + <^gqj,2 + <^gq,f,3 



II I II 
^gq,f,l + '^gqj,2 

JI 



a 



JI 



gq,error ' gq,error 



(5.11) 



—Cqg^s / dx (3{sx) I dy- / dcos9[ I d62 Ff^ (s,x,y, cos 6[, 62) 

' J p{s) J —1 1 — y J —1 Jo 

-Cggs dx (3(sx) C dy^— C dcos9[ f dO'^ Ff (s, x, cos^i, 
' Jp(s) Ji-yo i — y J-i Jo 

{x' + (1 - xf) 



lp{s) 

f 

lp{s) 



2 I A, 



dx 

y/ "Pis) X 



X dcos9[Y,\M'^^'"^^ixs,xb+)\^^^^ 

C,„3 A, j^^^^ dx P{sx) d cos 9, d92 



y=l 



+0(AJ) 



Pis) X 



+ {x' + (1 - xf) 



2 2 \syo, 



ln(l - x) 
d cos 9[ 



2(0) 



-Cqg,3 P{s) dz dy dv- / dcj)^ Ff^ {s, z^y^v^cj)^) 

Jo J-1 J-1 1 — V Jo 

-Cqg,3 (3{s) dz j^^dy dv- / rf02 Fff {s, z,y,v,(/)'2) 



vo 



-Cqg,2 epn i^-^J (3{s) dz 
Cgg^s A^ J^z J'^dyj^ 

l + {l-z)^\ 



1-A„ I rn 

dv / 

VQ 1 — V Jo 

n + {i-zf' 



z 



'II 



2(1) 



,/ f dFl]{s,z,y,v,(t)'2) 
dv 



+ O(A^) 



v=l 



Cqg,2 4 /3{S 



X J dz iz 



SVq 



z^{l- z)^p^{s) 



X 



j\^gq-*Wq 



S. 



'II 



2(1) 



drf,q{r) J DQ,Dq2^^-^^ [g(pi)q(P2) ^ W{Q,)q(Q2)] 

Cqg,4P{s)J^ dTf^qir) d COS ^1 ^ 



T 



(5.12) 
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The comments after (|5.4|) are also valid here. 

All the above terms will contribute in the 2 body inclusive scenario and in the jet 
scenario. However, in the 1 jet scenario, as we mentioned in section [IV. D| , by setting xq = 1 
and yo = vq = we are left only with the following contributions 



I II ~ p 

^gg = %,/,! + %,/,! + ^qg 

^OT = + ^OT- /,i + '^M • (5-13) 
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